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Abstract: The focal point of this paper is to provide some simple and efficient criteria to 
judge the I?-stabihty of two families of polynomials, i.e., an interval multilinear polynomial 
matrix family and a polytopic polynomial family. Taking advantage of the uncertain parameter 
information, we analyze these two classes of uncertain models and give some LMI conditions for 
the robust stability of the two families. Two examples illustrate the effectiveness of our results. 
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1 Introduction 


The study of robust stability problems under parameter uncertainties has been pioneered 
by the Russian scientist Kharitonovf 1978). A rich array of useful results have been developed 
over the last twenty years0“0’0“liiL Generally speaking, by dealing directly and effectively 
with the real parameter uncertainties in control systems, we can identify apriori the critical 
subset of the uncertain parameter set over which stability will be violated. The seminal theorem 
of Kharitonov\^ points out: any real parameter interval polynomial family is Hurwitz if and 
only if a special subset (called Kharitonov set) is Hurwitz. To general uncertain systems, edge 
theorem gives a positive answer®, which is an one-dimensional test. 


Consider the unity feedback system with an interval plant and a fixed controller in for¬ 
ward path, its characteristic polynomial is a multilinear function of certain interval variables 
0. That is to say, when all but one variables are fixed, all coefficients of the poly¬ 
nomial are affine linear in the remaining variable. The collection of all such models is called 
multilinear uncertainties structure. For an MIMO system, if all relationships between any input 
and any output belong to corresponding polytopes (the simplest form is a line), all admissible 
models form a polytopic polynomial matrix family. In the pas t one o r two decades, there is 
continually growing interest in the robust analysis of matrix® H’ 11®’ 1111’ 11® . Unlike polynomial 
case, the vertex result does not hold for interval matrix. In fact, there does not even exist a 
result similar to Edge Theorem for general matrix. Many problems sti ll remain open until now. 
Work reported to date shows that a reduced dimensional test holdsll® for polytopic polynomial 
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matrices. Limited by the complexity of such problems, several methods have been proposed, 
such as eigenvalues estimation, Lyapunov approach, algebraic approach and spectrum theory, 
etc. §-i. Usually, the algebraic approach based on Kharitonov theorem is hardly effective 
and convenient when used for the robust analysis of matrix directly. Lyapunov approach is an 
appealing method developed in the context of Lyapunov theory, which presents in the form of 
Linear Matrix Inequalities(LMI). However, until now there is few useful result on robust test for 
matrix families under parameter uncertainties. 

The purpose of this paper is to address the P-stability of interval multilinear polynomial 
matrix family and polytopic polynomial family. Even though they are nonlinear problems, we 
still can establish some efficient robust stability tests, which are usually negative for a general 
nonlinear case. Rec^t work addressed the T>-stability problem for polytope of matrices using 
Lyapunov approach® and spectrum theory®. In this paper, we adopt both of two methods to 
analyze two classes of uncertain MIMO models and give several LMI criteria on robust stability. 
Our technique not only captures the uncertain parameter information, but also is easy to use. 
In the end two examples demonstrate the effectiveness of our results. 

2 Definitions and Notations 

In this paper, we use the following standard notations and definitions. 

Definition 1 Given an open convex region D in the complex plane, a scalar matrix is termed 
D-stable, if all its eigenvalues lie in D] and a polynomial matrix is termed D-stable, if all roots 
of its determinant lie in D; a matrix family is termed D-stable, if all its members are D-stable. 

Definition 2 Let M be an arbitrary set, we define convM as the convex hull of M, i.e., the 
smallest convex set which contains M. 

Definition 3 For a matrix A, its right null-space Ma is defined as the space whose every ele¬ 
ment Na satisfies ANa = 0. For simplicity, we denote A/a a basis for the right null-space of A. 

Definition 4 A polynomial matrix is a matrix with all of its entries being polynomials; an inter¬ 
val multilinear polynomial matrix is a polynomial matrix with all of its entries being multilinear 
dependent on some interval coefficients; a family of such matrices is called interval multilinear 
polynomial matrix family, such as the model in (5); a polytopic polynomial matrix is a matrix 
with all of its entries being polytopic polynomials; a family of such matrices is called polytopic 
polynomial matrix family, such as the model in (6). 

Definition 5 Let P C C be an open convex set of the form 



where H is an 2 x 2 matrix and B* = B. V is called an LMI region(®]||T^). 

Definition 6 For every i G A^}, Aj(s) is an n x n polynomial matrix of the form 


Aj(s) — Aq -|- A’^s a 



( 1 ) 






and the n x nl scalar matrix Ai = {Aq, ..., Aj) is the coefficient matrix of 

Definition 7 For every i G {1,..., n}, j G {1,... , n}, Vij{s) is a polytope of polynomials , i.e., 

{ m m 'j 

J2XkP^\s): Afc>0, ^Afc = l , (2) 

fc=i fc=i J 

where p\^\s),... ,p\'p\s) are fixed /cth-order polynomials. Apparently, the vertex set of every 
Vij{s) is 

Kiji-s) = {p[f{s), k = l,...,m} (3) 

and the exposed edge set of every Vij{s) is included in the following set 


E: 




(s) = + (1 - X)pf^{s), k,t = I,- ■ ■ ,m, k^t, A G [0,1]| 


(4) 


Definitions Denote Q = |(gi,: O'* G ai{qi,... ,qm), ■ ■ ■ ,aN{qi, ■ ■ ■ ,qm) 

are multilinear functions of qi,... ,qm- 


Definition 9 


( I 


R = 


O \ 


I O 


O I 


\ o I J 

is a 2nl x (1 + l)n-dimensional scalar matrix and I, O are nxn unity, zero matrices respectively. 
As usual, (8) stands for the Kronecker product. 


3 Preliminary results 


In what follows, two kinds of uncertain models are considered: 

MA{s) = |^ai(gi,... ,gm)k('S), (o'!, • • •,G q| 

VA{s) = {{pij{s)) 

nxn : Pij{s) G Vij{s), i,j = 1,... ,n} 

The vertex sets of XAA{s) and VA{s)^ respectively, are 

Kma{s) = {Eiliaj(gi,...,gm)k(s), qj G {q^,q^}A = 

Kva{s) = {(Pij('S)) nxn : pij{s) G Kij{s), i,j = 1,.. . ,n}. 

Let P” be the collection of all permutations of 1, 2,... , n, and define 

p ^ Pki^i-s) e Eki^{s), e P;^, k = l,...,n 

^VA[S) - < [mAs))nxn . eKk^,is), ik = I,...,Ik-1,Ik+ l,...,n 


(5) 

( 6 ) 

(7) 

( 8 ) 


It is easy to see that, Ep_a{s) is a subset of VA{s) produced by taking only one entry from its 
exposed edge set in every row/column and all other entries from their vertex sets. 



The lemma below is due to Henrion, et al. 

Lemma iQ A polynomial matrix Aj(s) is stable if there exists a matrix Pi solving the LMI 
feasibility problem 

MXW{B ® Pi)RMA, < 0 , P^ = P*> 0 . 

Where is the right null-space of Ai, and * denotes the conjugate transpose operator. 
Lemma 2 ||^ (Edge Theorem) Suppose P C C is a simply-connected region, for any polynomial 
polytope n without degree dropping, the root set of O is contained in P if and only if the root 
set of all exposed edges of Q is contained in V. 


Another lemma is on the P-stability of the family VA{s), 

Lemma 3 |p^ PA{s) is P-stable if and only if E'pj\^{s) is P-stable. 

Proof: For all € PA(s), let ^*^( 5 ) = {Pijis))nxn^ where Pij{s) G P{s)- For simplicity, we 

write pij for every Pij{s) G Vij{s). In what follows, we will construct several sets in terms of 
A°(s). Let Ak = {Afc(ii ,... ,in\ s), ii,..., 4 G { 1 , ... ,n}} (/c = 0,...,n), where 
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It is easy to see that A°(s) = ^0 and Ak C Ak+i- In the sequel, we will prove our statement in 
two steps: 

1) Firstly, we will show that An is P-stable if Epj^{s) is P-stable. By definition, we have 


f • • • ? ini '^) 

qit G 1 

-^itt ('5)) i — it 
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For all An{ii,... ,in] s) G An, if {ii,...,in) G -P", then A„(ii,..., s) G P-p^(s). Otherwise, 
there must exist some pair is,it satisfying is = it. Without loss of generality, suppose ^1 = ^2 = !, 
namely 


/ 


An (il, ■ ■ ■ ,in, s) — 


^11 912 


91^ 


9ii G Pii(s) 912 G £'12(5) 


y 9nl 9n2 ••• Qnn 

By using Laplace formula on the first row of An{ii,... ,in', s), we have 


n 

det{Aniii, ...,in, s)) = 9ii^ii + 912 M 12 ^ quMu 

i=3 


where Mu is the algebraic complement of qu. By Lemma 2, 

A„(zi,... s) is P-stable g'nMn-h 912^12 + EL 3 9 ii^ii and 

qiiMu + 912 M 12 -h T,i =3 quMii are P-stable. 










The corresponding matrices are 


A 

gi2 • 

- gin ^ 

and 

^ gii 


■ gin \ 

Inl 

qn 2 ■ 

• Qnn y 


^ gm 

qn 2 - 

qnn y 


where € Ku, (?i 2 G ^ 12 ■ For these two uncertain matrices, if they do not belong to Epa{s), 
then there must exist at least two equal indexes. Repeat the same process to them, in the end, 
we have 

E-pa{s) is P-stable => An is D-stable. 

2) Secondly, for all An{ii■,■■■■, in'-, s) G An-, by using Laplace formula on the n-th column of 
An{ii-, - - - -,in'-, s) and by Lemma 2, we have 

An is P-stable => An-i is T>-stable. 

Continuing this process, we have 

Ak is P-stable ^ Ak-i is T>-stable. 

Since 1) and 2) imply that is D-stable. That is to say, if E-p^is) is P-stable, 

then, for all G VAis), is P-stable. By definition, this means that VA(s) is D-stable. 

This completes the proof of Sufficiency. 

Necessity is obvious because E-p^is) is a subset of VA{-s)- o 

In this paper, we assume that both VA{-s) and AiA{s) have fixed degrees. 


4 Main results 

Although overbounding is a bit conservative, it still offers a powerful tool to solve the stability 
problem for control systems with multilinear uncertainties. 

Theorem 1 M.A{s) C conv{KA 4 A(s)}- 

Proof: For any Ao(s) G M.A{s)., by induction, we will show Ao(s) G conv{K_MA(s)}- Denote 
q = {qi, - - - -,qt, • • •, <?m)^ G Q, where qi,... ,qt are interval parameters and qt+i, - - - ,q-m are 

fixed. 


If t = 1, i.e., qi G [gf, c^] and q 2 , - - - ,qm. are are fixed. Then, in this case, 

MA{s) = [j2f=iai{qi,---,qm.)Ai{s), qi G 
Kma(s) = [j2fLiai(.qi,---,qni)Ai{s), qi G {qi,qi }^ - 

Since ai(gi,..., g^), • • •, aAf(gi, • • •, gm) are linear in gi, ai(gi, • • •, qm)Ai{s) is also linear 
in gi. Clearly, Ao(s) G conv{KA4A{s)}- 


Assume that the claim holds for t = k. When f = A: + 1, we have, in this case, 

qi e [qi,qY]A = i,...,A: +1 1 


MA{s) = <! Eili a*(gi, • • •, gm)Aj(s), 


gfc+ 2 ,..., gm are fixed. 


X/i=l (i'iiqii ■ ■ ■ 1 qra)AiiyS) 1 


qi e {g» ,g» }d = + i 

gfc+ 2 , -- - ,qm are fixed. 



For all Ao(s) G M.A{s), there exists an m-dimensional vector := , q%)'^ G Q, satisfying 


N 


^o(s) = • • • ’ 


2=1 


Since ai(gi,..., g^), • • •, am(gi, • • •, <?m) are linear in qu+i, we have 

aMi, Qk+i’ ^oa*(g?, • • •, g^+i, g°+2> • • • > 9m) 

+(1 - Ao)ai(g?, ...,ql, •••,9m) 

for some Aq G [0,1]. Therefore, 

^o(s) = E)Ii^(Aoai(g?,..., gfc+i, gfc+2,..., O + (1 - Ao)a*(g5’,..., g^+i, 9^+2, • • •, 9™)) ^*(s) 

= Ao E)Ii ai(gi, • • •, 9°, 9fe+i, 9fc+2, • • •, 9m)^i(s) 

+(i - Ao) E)Ii ai(gi, • • ■, 9fe, gl+i, gfc+ 2 , • • •, q^)A{s) 

Since both X;ili ai(9i, • • •, 9fc, 9fc+i, 9fc+2, • • •, 9m)^i(s) and «i(9i, • • •, 9°, 9^+i, 9^+2, • • •, 9m)^i(s) 
belong to the set 

f /- ^ /I ^ 9i £ { 9 /", 9 F },^ = 1,---,^ 1 

9 fc+i, ••• , 9 m are fixed. J 

we conclude 

/!„(»)€ ^ (rf • iT. • = 1. , t + 1 1 

9 fc+ 2 ,..., 9 m are fixed. J 

That is to say, the claim holds also for t = k + 1. Therefore, our conclusion is verified inductively. 


For all Bi{.s) G K_\ 4 j({.s), rearrange it as Bi{.s) = Bq + B\s+. . .+BisK Take Bi = {B^, ..., Bf) 
the coefficient matrix of Bi{s). It is easy to see that there exist 2™ distinct Bi{s). Hence, 

conu{iir_ 4 (s)} = conu{Hi(s),... ,H 2 m(s)}. 


With Lemma 1 and Theorem 1, we get an LMI condition for robust stability of interval multi¬ 
linear polynomial matrix family A4A{s): 

Theorem 2 A4A{s) is robust P-stable if there exist some matrices Pi = P* > 0,Q solving 
the LMI feasibility problem 


R ' 

* 

B(E) Pi 

Q ' 

’ R ' 

Bi _ 


Q* 

0 

. . 


i = 1,...,2™. 


(9) 


Proof: For every H(s) G MA{s), by virtue of theorem 1, there exist Ai,..., A 2 "* G [0,1] such 
that J2f=i Aj = 1 and H(s) = Ym=i KBi{s). Moreover, for all i G {1,..., 2™}, 


R ' 

* 

B®Pi 

Q ' 

' R ' 

Bi _ 


Q* 

0 

. . 


^R*{B(^ Pi) R + B*Q*R + R*QBi < 0 

^ El=i Ai {R*{B 0 Pi)R + B*Q*R + R*QBi) < 0 

^R*[B(^ (Eti ^^Pi)) R + (Eti KBi)*Q*R + R*Q{Ylli hBi) < 0 



Multiplying from the right and from the left, the inequality becomes 

R 


N" 


A 


R 

* - 




Q* 0 


EI=1 


A/U < 0 


Because of M_a = A/’^ 2 ™ we have 

Z^i=l 


/ 2 ’" \ 

MXR* A*^*) BMa < 0 


i=l 


From P* = Pi > 0, we have that Aj-Pi = (Z]i=i AjP^)* > 0. Now by the Lemma 1, the 
conclusion is obvious. o 

Remark 1 The standpoint of Theorem 2 is to transform stability problem into a positive 
real-like condition, and the latter can be solved using the LMI toolbox. 


In regard to Lemma 1 and Lemma 3, we claim that the stability of uncertain family PA{s) 
can be inferred from whether an LMI condition holds or not for K-pa{s)- This is shown by the 
following two theorems. 

Theorem 3 VA{s) is D-stable conv{K-pA{s)) is I^-stable. 

Proof; Sufficiency; We will show E'pa{s) C conv{K'pA{s)). 


For all A-i(s) G £'-p_ 4 (s), by the definition of £'-p^(s), there exists (Zi ,... ,ln) & Pn such that 
A .(„\\ . Pklki^) ^ -P'fcZfc(s), Zc = 1 ,..., n, 

^IV^/ — ) Jnxn • / \ ^ \ ’ 1 7 i7ii 

Pkik (^) ^ ('^j Ik — 1; • • • ; 1, • • • , 

For pi/j (s) G Eii^ (s), we know that there exist m real numbers An,..., Aim & [0,1] satisfying 


Pihis) = AifcPui Aifc = 1- 

k=l k=l 


Using addition of matrices, we have 


Ai(s) — ^ XuA 


(t) 


t=i 


where is the matrix that all its entries coincide with vli(s) except one, which lies in the 
first row and the Zi-th column and equals to for every t G {1,..., m}. Thus, for everyone of 
..., all of its entries in the first row belong to vertex sets. 


For every A^^^ {t = 1,... ,m), applying the same process to P 212 , we can find m uncertain 
matrix families, and for every matrix which belongs to one of those families, all its entries in 
the first and second rows belong to the corresponding vertex sets. Continuing this procedure, 
we will get Ai(s) G conv{K-pA{s))- Therefore, P-pyi(s) C conv{K'pA{s)). Then, 


conv{K'pA{s)) is ID-stable Ep^is) is T’-stable 
( by Lemma 3) VA{s) is T>-stable. 



Necessity: The relationship between 'PA{s) and conv{K'pji,{s)) can be easily established, 
thereby Necessity is proved. For all ^i(s) G conv{K-p_A{s)), there exist numbers Xij G [0,1] 
J27,j=i 1 matrices Fij{s) G such that 

n 

^l(s) = ^ijF'ijis)- 
hi=i 

Denote Fij{s) = (fhAs)) with /AAs) G Khi{s) for all h,l G n}. By addition of 

V / nxn 

matrices, 

^i(s) = ( £ 

Vhi—1 / nxn 

For every h £ {1, • • • ,n},l G {1, • • • ,n}, YAi,j=i ^ijfhii^) still belongs to Vhi{s) whenever Xij G 
[0)1]) YAi,j=i^'i -3 ~ 1 ^ Fhi{s). Thus conv{K'pj[{s)) C VA{s). This completes the 

proof. o 


For all ^(s) G VA{s), we have T(s) = {pij{s))nxn, where Pij{s) G 'Pij{s) with degree 1. 
Rewriting it as 

T(s) = Aq + Ais + ... + AisK 

Denote A = (Aq) • • •) A)) then ^ is an n x n/ scalar matrix. By Theorem 3, we have 
Theorem 4 VA{s) is robust D-stable if there exist some matrices = PJ > 0, Q solving the 
LMI feasibility problem 


’ R ' 

* 

B®Pa 

Q ' 

’ R 

A 


Q* 

0 

A 


< 0 , A£K'p^{s). 


( 10 ) 


Proof; By Theorem 3, this problem is equivalent to the P-stability of conv{K'p_A^{s)). For the 
latter, for all vl(s) G conv{K-p_A{s)), there exist n? numbers Xij G [0,1], — 1 ^'^‘1 

matrices Fij{s) G P-p^(s) such that 


T(s) = ^ A,,P,,(s). 

hi=i 

Now by a similar argument as in the proof of Theorem 2, we get the desired result. o 


5 Illustrative Examples 

In this section, we give examples to illustrate the utility of our main results. Example 1 is 
considered in the context of robust stability of interval multilinear polynomials with respect to 
left half plane. 

Example 1 (n=l, N=2, 1=3, m=3) Let ^i(s),yl 2 (s) be two given polynomials 

Ai(s) =s3 + 2.64s 2 + 1.82s+ 0.37 (11) 

^ 2 ( 5 ) = s^ + 5.57 s 2 + 9.04s+ 3.85 (12) 

And the uncertain model is Al(s) = {ai(gi, § 2 , (? 3 )Ai(s) + 02 (( 71 , 02 , 93 )^ 2 (s)}, where (71 G [1,2], (72 G 
[3,3.8], qa G [0.5, 0 . 8 ] and 01 (^ 1 , ^ 2 , 93 ) = O.Og'i + 0 . 1^2 - g's + 0 . 191 ^ 2 ) 02 ( 91 , 92 )®) = - 0 . 6 ^ 1 - 



0 . 1 g 2 + gs + 1 — 0 . 01 g 2 <? 3 - For Hurwitz stable, the 2x2 Hermite matrix B corresponds to 

^ Q . Applying Theorem 2 to this problem, it suffices to solve 16 linear matrix inequalities. 

Using the LMI Toolbox in Matlab, then it is easy to check that the corresponding LMI problem 
is feasible. Thus, we conclude that the whole polynomials family is robust Hurwitz stable. 
Example 2 (n=3, N=2, 1=3, m=3) Consider the third order uncertain model below 


A(s) = {ai(gi, 92,93)^1(5) + 02(91,02,93)^2(5)} 
^i('S) = 

/ 15s 3 + 2.55^ + 12s - 1 
10s2 - 7.4s + 2.8 
12 s 2 - 1.1s -4.7 


10 


V 

^ 2 ( 5 ) = 

( 20s3 - 13s 2 - 18s - 0.96 
-9.8 s 2 + 14s - 4.5 
y 6s^ + 17s + 7.1 


13s3 + 12s 2 +0.22s-1.1 
23s3 + 4.7s2-1.3s + 3 
2.8s 2 + 0.99s - 0.5 


2.7s3 - 7.3s 2 - 7.9s - 3.7 
11s3-16s2 + 15s + 7.6 
25s3 - 2.2s 2 - 3.8s + 0.035 


-3.7s3 + 13s 2 + 0.078s - 6.7 
19 s3 + 0.35s 2 - 7.5s + 3.4 
-1.3 s 2 + 1.1s + 1.5 


-20s3-4.5s 2 +9.4s-0.27 
-11s3-5.5s2-20s + 3.8 
lOs^ - 1.9s2 - 10s - 5.1 


with 9 i G [1,1.2], 92 G [2.1,2.4], 93 G [1.5,1.8] and 01 ( 91 , 92 , 93 ) = 9i - 92 + 93 + O.I 9192 , 
02 ( 91 , 92 , 93 ) = —91 + 92 — 93 + 1 — O.OI 9293 . In this Example, the quadratic stability region is 


the unity circle, thus the associated matrix is 


-1 0 

0 1 


Solving the corresponding LMI in 


Theorem 2, it is easy to show that the uncertain model is robust Schur stable. 

Remark 2 Our results can also be verified by the plots of root loci of the whole polynomials 
family in following figures. From the plots of root loci, we can see that our LMI criteria are not 
very conservative, and can provide correct, effective information on robust stability of uncertain 
systems. 


6 Conclusion 

In this paper, we have dealt with the performance robustness of interval multilinear polyno¬ 
mial matrix families and polytopic polynomial matrix families. Some computationally tractable 
and nonconservative sufficient conditions for these two classes of system models have been ob¬ 
tained. 
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